Short wavelength spectrum and Hamiltonian stabihty of vortex rings 
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We compare dynamical and energetical stability criteria for vortex rings. It is argued that vortex 
rings will be intrinsically unstable against perturbations with short wavelengths below a critical 
wavelength, because the canonical vortex Hamiltonian is unbounded from below for these modes. 
To explicitly demonstrate this behaviour, we derive the oscillation spectrum of vortex rings in incom- 
pressible, inviscid fluids, within a geometrical cutoff procedure for the core. The spectrum develops 
an anomalous branch of negative group velocity, and approaches the zero of energy for wavelengths 
which are about six times the core diameter. We show the consequences of this dispersion relation 
for the thermodynamics of vortex rings in superfluid *H[e at low temperatures. 



I. INTRODUCTION 

The long wavelength oscillation spectrum of large vor- 
tex rings in incompressible, inviscid fluids is established 
since the pioneering work of W. Thomson (Lord Kelvin) 
§, J. J. Thomson @, and Pocklington |]. The vaHd- 
ity of that spectrum is restricted to wavenumbers much 
less than the inverse core size and rings which are large 
compared to the extension of the core. There are, how- 
ever, processes for which it is desirable to know the large 
wavenumber properties of the spectrum for smaller rings: 
Vortex ring nucleation, reconnection of vortex filaments 
and dissipation in the turbulent energy cascade are be- 
lieved to occur on very small length scales, reaching down 
to a few times the vortex core size. Because a line vortex 
represents a string object, elasticity modes will be ex- 
cited during the rapid movements executed by the string 
on small length scales. A fluctuating line should have 
equilibrium states different from a non-fluctuating one, 
because the quantum or classical statistical fluctuations 
renormalize the total free energy as compared to the un- 
deformed ring. 

In what follows, we shall derive the collective, small 
amplitude oscillation modes of a vortex ring in an incom- 
pressible, inviscid fluid. The dispersion relation is exact 
within the geometrical cutoff procedure we employ, and 
displays a maximum and an anomalous branch of nega- 
tive group velocity. The critical wavelength for the spec- 
trum to possess a positive excitation energy corresponds 
to one oscillation of line within a length about an order 
of magnitude above a geometrically defined core size. It 
will be argued that, due to the properties of this spec- 
trum and the structure of the canonical Hamiltonian, a 
vortex ring is potentially unstable in an intrinsic man- 
ner, because the Hamiltonian is unbounded from below 
for short wavelengths. Physically, the energetical insta- 
bility is caused by the fact that quantities playing the 
role of "mass" and "spring constant" in the Hamiltonian 
do simultaneously assume negative values. The energeti- 
cal instability, taking place for large perturbations which 



are of wavelengths less than the critical wavelength (of 
about six times the core diameter in our core model), 
occurs though the ring is dynamically stable for nearly 
all wavelengths down to the core size. The dynamical 
instability of a vortex ring, related to the ocurrence of 
imaginary excitation frequencies, is only relevant for cer- 
tain critical wavelengths. We will thus present in this 
paper an argument that the relevant stability criterion 
for a vortex should be that of energetical stability. 

Below, in the section to follow, we will first introduce 
an action principle which gives a transparent representa- 
tion of the vortex dynamical behaviour in incompressible, 
inviscid fluids, from which the canonical Hamiltonian for 



small perturbations in section III, representing the vor- 
tex eigenmodes, naturally follows. Section ^ gives an 
account of the thermodynamics of vortex rings, related 
to the excitation spectrum on a quantized vortex in su- 
perfluid "^He, where the consequences of the predictions 
made in this work should be particularly clearly seen. We 
conclude with some remarks. 



II. DERIVATION OF THE OSCILLATION MODES 
A. Action Principle 

The peculiarity of the dynamical behaviour of vortices 
in the incompressibility approximation, the fluid having 
a constant mass density pa, consists in the fact that con- 
flguration space and phase space coincide . The mo- 
menta are, in this limit, generally expressible as functions 
of the co-ordinates, and play no independent dynamical 
role. This fact gives rise to the following action functional 
of the line conflguration C, in terms of the positions R 
of line elements dR with a constant velocity circulation 

r i: 



S[C] = dt (^-^ Po ^ {dRAR, dfR) - H [C] 



(1) 
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The factor ^ in the kinematical term reflects our choice of 
cartesian co-ordinates in what follows, and corresponds 
to a (co-ordinate) gauge for the vortex momentum 0. 
The vortex kinetic energy is given by the Biot-Savart 
expression 



„r^n r2 / / 1 {dR,dR!) 



(2) 



where the shorthand notation R — R{(f>, t) and R! — 
R{(j)',t) is used. This relation for the energy yields the 
usual asymptotic logarithmic dependence of the station- 
ary vortex energy on the infrared cutoff L (system size, 
distance to the next vortex, or radius of a vortex ring) 
and the ultraviolet cutoff (the core size), in the form 
ln[i/(C^c)], where C is a core model dependent constant. 
Stationarity of the action for first order variation of the 
action (0) after R leads to the local velocity of a line ele- 
ment being perpendicular to the line element, and given 
by the Biot-Savart nonlocal induction law: 



dRA\ dtR' 
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The above proves that the action (|I|) leads to the correct 
equations of motion, familiar from the fluid dynamics 
literature 

Let R {(j), t), for < </> < 2tt, describe the instantaneous 
shape of a moving vortex ring at time t, fluctuating with 
amplitude u((j), t) around its circular equilibrium shape 
Ro{(t),t), so that R{(t>,t) = Ro{(t),t) + u{(j),t). Then, 
we parametrize line element position, equilibrium posi- 
tion, and small perturbations around this equilibrium as 
follows 

R{(j),t) ~ R±{(f>,t) {cos(f>ex + sin(f>ey) + R\\{(f>,t)ez 
R± =ro + u±{(f),t) , i?^ II = + w II ('/', • (3) 

From this choice of co-ordinates and the form of the ac- 
tion (|l|), the phase space variables for small oscillations 
tt(0, t) are concluded to be 



qi(t>,t) = Mil {(j),t) 
p{(f>,t) = {Tporo)uj_ {(t),t) 



(4) 



These phase space variables are employed for the descrip- 
tion of the vortex eigenmodes which follows. 



B. The spectrum 

The fundamental cutoff to be introduced for the con- 
tinuum description to be valid is that the separation of 
two line elements should always exceed a length ^c, 



\R-R'\ > 



(5) 



The length is thus defined as the cutoff diameter of 
the vortex core. The above prescription is the simplest 
exact procedure to ensure that the Biot-Savart integrals 
remain regular. If the Biot-Savart description is refined 
by, e.g., a density profile in the core smoothly increasing 
within a distance to the constant po, instead of being 
cut off to be exactly zero at this will effectively 

yield a different ultraviolet cutoff, that is, a different core 
constant C of order unity, multiplying a (fixed) value 
of ^c- However, the dynamical behaviour of the vortex 
line on a scale well outside the core domain will not be 
affected by the core model. 

To leading order in the fluctuations, the condition (||) 
is equivalent to 



> 



2ro 



(6) 



Introducing the above condition in the integrals deter- 
mining the velocity in (|^) by means of a Heaviside step 
function, we obtain the equilibrium velocity of the ring 



Vo = 
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In [cot ((5/4)] , 



(7) 



where the cutoff angle is determined by the parameter 



S — 2 arcsin . 

2ro 



(8) 



To obtain the ring oscillation modes, we expand the small 
quantities U|| and u±_ in a Fourier series, u±{4>,t) = 

E„"-L.nW e'"*, ui\i^,t) = E„W||,„(t) e™"^, and use 
the above described cutoff procedure of ^ respectively 
(^). We then obtain the linearized equations of motion 
for parallel and perpendicular oscillations of the filament 



The coefficients in this linearized version of 

r 1 
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r? In cot 1\ 



(9) 



(10) 



1 + cos(n(5) , 2\ 1 T 
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are given in terms of integrals containing, due to our 
parametrization (|) of the ring oscillations, trigonometric 
functions only: 
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In Fig. 0, the coefficients a„ and 6„, for tq/^c = 7.5, 
are shown. 
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FIG. 1. The coefficients a„ (stars) and fen (boxes) in the 
equations of motion in units of loc defined in Eq. (p^, 
for ro/^c = 7.5, up to n = 12. There occurs a dynamically 
unstable mode of imaginary frequency for n = 8. 

According to the frequencies of oscillation of a vortex 
ring evaluate from 



(12) 



We scale the coefficients a„, bn in Figs. and frequen- 
cies LOn in Fig. ^ below in units of the fundamental cy- 
clotron frequency of the vortex core. 



(13) 



the frequency with which a point on the core revolves 
around the line designated by i?(0, t) (the name stem- 
ming in the magnetic analogy from the role of F as a 
flux). 

The frequencies (^2|) are exactly zero for both n ~ 
and n = 1, to any order in 5. In the first case, oq = and 
6o = (r/47rrQ)(cos((5/2))~-^ — vq/tq, whereas in the latter 
case, ai = vq/vq and hi = 0. The first case of symmetry 
corresponds to u±^o =const., and tells us that the radius 



of the ring as a function of i5 is determined up to the 
(constant) value of u±fi. In the second case, in turn, uy ^ 
is a constant. The resulting line deformation resembles, 
for (5 — > 0, a translation of the ring as a whole. Except 
for radii tq which are just about an order of magnitude 
above ^c, the coefficients a„ and 5„ are practically the 
same over the range of allowed values of n, so that for 
i5 <C 1 the frequency ujn is essentially equal to either 
|a„| or |5„|, and the waves around the ring are circularly 
polarized, like for ordinary Kelvin waves. However, for 
To getting closer to ^c, dn becomes increasingly different 
from bn, and the waves become elliptically polarized, the 
absolute ratio of amplitudes in the ring plane and out 
of the plane being given by |uj.,„/u||,„| = y^\a„/bn\. A 
small ring thus oscillates more in the ring plane than 
out of the ring plane. Let us also stress that for the 
calculation of the oscillation modes, the nonlocality in 
the Biot-Savart integrals of (||) respectively is fully 
taken into account. 




FIG. 2. The coefficients a„ given in Eq. (^o|), up to 
n = 2tt5~^ ~ 27rro/^c, for ro/^c = 60, in units of lOc de- 
fined in Eq. (^|). For these small values of 6, a„ and b„ are 
indistinguishable within the Figure's resolution. The dotted 
curve is the asymptotic result for 5 ^ in equation ([l4|). 

For direct comparison with the dispersion of Kelvin waves 
on rings, we consider, for fixed n, the limes (5 ^ in ( p^ ) 
and (pT|), to obtain 
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(14) 



where 5„ — ^ 1)^^- There is an important dif- 

ference of the dispersion relation (^2|) (which is exact 
within our hollow core model), with a„ and bn from ([l^), 
to the usually quoted asymptotic results of Lord Kelvin, 
J.J. Thomson 0J| and Grant (also cf. the work of 
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Pismen and Nepomnyashchy , which found the same 
result as Grant, but within a much simpler scheme, simi- 
lar to ours) . These last results correspond to the relations 
( p^ ) for the coefficients a„ and 6„ in the limit ofS ^ for 
fixed n (save for different core structure constants). The 
important difference consists in the fact that the geomet- 
ric cutoff prescription (j^), which ensures that a core of 
diameter is always excluded in the evaluation of (^, is 
taken care of in relations ( p^ ) exactly for any admissible 
value of the ratio of ring radius and core diameter ?'o/Ccj 
such that n5 ~ 0(1) can be consistently realized. The 
anomalous branch also occurs, shifted to smaller mode 
numbers, as a consequence of relations (14). However, 
the minimum resides at values of n which are 

beyond the applicability of equations (^. We have de- 



corresponds to a transformation to the rest frame of a 
ring of radius tq, moving with velocity vq. The above 
expression for Ti, then, represents the rest frame Hamil- 
tonian of the vortex. The phase space variables may, for 
example, be chosen to be g„ = wy „ and Pn = Po^'i'ou±,n, 
so that the mass M„ = poTro/bn and elasticity (spring) 
constant D„ = poTroCLn- Equally well, we may choose the 
option g„ = M_L,n and p„ = —poTrQU\\,ni which reverses 
the role of mass and elasticity coefficients in conventional 
Hamiltonian language (replaces a„ by 6„ and 6„ by a„ in 
M„ and Dn). The identity of phase space and configura- 
tion space P R| implies that both options are viable. 



picted the difference between the exact and asymptotic 
results in Fig. ||, for the whole range of n up to the 
value n = 2^6'^ ~ 27rro/^c, corresponding to /c ~ 27r/(^c 
(A ~ ^c)- For the value tq/^c — 60 used, the coefficient 
bn is essentially identical to a„ within the resolution of 
the Figure (cf. Fig. |], which has tq/^c — 7.5, and where 
the difference between a„ and bn is clearly discernible). 

For completeness, we state the solution of the equa- 
tions of motion (^). If the vortex ring undergoes at time 
t — a, deformation represented by 



zill (0,0) =5R5]4„e™'^ 

n 

Ui(0,O) =^J2 



(15) 



the solution at a later time t takes the form 

M|| (0, t) = ^ U|| „ COS (uJnt) + u° ^„ — sin (tj„t) 



u± I 



, t) = ^ cos {uJnt) - Uji_„ — sin (tj„t) 



Jn<t> 



(16) 



where ujn 
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III. HAMILTONIAN 



The Hamiltonian corresponding to the canonical vari- 
ables (||) and the action (|^) assumes the form 



Po^ro 



where the stationary energy of the ring is given by 



■(■■■ 






'S' 




(i). 


— 2 cos 




) 








2 





(17) 



(18) 



The (constant) variables wj_,o and um do not appear in 
( p7| ) because of our choice of parametrization (^, which 
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FIG. 3. Stable oscillation frequencies of a vortex ring, in 
units of the cyclotron frequency coc, as a function of the mode 
number n, for the ratio ro/^c = 60 ~ , up to Uc — n/{35). 
For this ratio of radius and core diameter, the exact value is 
ric = 63. The maximum is, essentially independent from the 
value of S, situated at max[ci;„] ~ 0.011 cJc, with n ~ nc/2. 



From the Hamiltonian (|lj), we gather that stable oscilla- 
tion modes are those which have a„ and 6„ both positive. 
They contribute positive energy to the Hamiltonian. En- 
ergetically unstable, though giving a real frequency, are 
modes which have a„ and bn both negative, because they 
contribute negative energy to the Hamiltonian. A dif- 
ferent sign of a„ and 6„ leads to dynamically unstable 
modes, which have imaginary frequencies and amplitudes 
u exponentially growing (or decaying) in time. For mode 
numbers above 



TT 

36 



(19) 



and up to n ~ 5S~^, both a„ and 6„ become nega- 
tive, such that the energy contribution corresponding to 
these modes is negative. To quadratic order in the oscil- 
lation amplitude, oscillations with wavelengths smaller 
than A ~ 6^c thus imply that vortex modes of such small 
wavenumber are unstable. Hence, the stable spectrum 
is restricted to mode numbers of magnitude less than 
ric — by definition the last mode number for which 
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the oscillation energy is positive semi-definite, before en- 
tering the negative energy domain seen in Fig. |^. We 
have plotted the dispersion relation of the stable modes 
in Fig. I for ro/Cc = 60. 

With regard to the validity of the assumption of an 
incompressible fluid, we note that the frequency at the 
maximum in Fig. situated at n ~ nc/2 — {2S)^^ 
for all values of ro/<^c not too close to unity, scales as 
ujn — 0.011 LJc, with the cyclotron frequency uJc defined 
in (|l|) (wc lO^^ggc-i in superfluid ^He (helium II)). 
Oscillation velocities thus remain, for moderate oscilla- 
tion amplitudes of the order of a few ^c, well below the 
speed of sound even in the superfluid helium II, where 
is of atomic size and CsS,c ~ F = i:lJc{£,c/'^)'^ , so 
that the incompressibility approximation holds. This is 
more questionable for the second, much larger, frequency 
maximum at n ~ 3<5~^, corresponding to the maximum 
negative value of the coefficient a„ in Fig. |[ given by 
a„ ~ — O.lScjc. 




FIG. 4. Shape of helical wave traveling along the ring, for 
mode numbers near the critical mode number ric, showing the 
vortex core displacement and the volume exclusion caused by 
this displacement. 

Consider, for a physical interpretation of the energetical 
instability around n — ric, Fig. ^ where we show the 
shape of deformation of the vortex core for a small wave- 
length of order Ac — 6£,c, corresponding to the crossover 
to the unstable oscillations regime. We may infer that the 
negative oscillation energy, occuring at a smaller wave- 
length than that shown in Fig. ^ is due to a volume 
exclusion effect. The excluded core volume kinetic en- 
ergy is large enough such that there is only small en- 



ergetical cost of exciting a perturbation on the filament 
for the stable modes with mode numbers slightly below 
Uc, and an energetical gain for the unstable ones. In the 
Hamiltonian (0), we neglect the tail of positive excita- 
tion energies corresponding to positive values of a„,6„ 
seen in Fig. ^, because it is very close to the limit of 
core elements touching themselves, at which point our 
core model certainly becomes invalid, because it is then 
meaningless to speak about oscillations of a hollow torus. 
That the coefficients a„ and &„, and thus the excitation 
energy, do increase again after n ~ 35~^, can be traced 
back to the fact that line elements having like circulation 
approach each other closely if we further compress the 
spiral of Fig. [i] along its axis. The volume energy ex- 
clusion effect we just described is then counterbalanced, 
for these very short wavelengths, by the resulting strong 
repulsion of adjacent elements of the same circulation. 

The frequency ( |l2[ ) is imaginary, if a„ and &„ have dif- 
ferent signs, and a dynamical instability results |]ll| , |l2| . 
We stress, however, that the unboundedness of the 
Hamiltonian ( [l7| ) from below leads to the energetical in- 
stability of the ring for mode numbers beyond Uc- This 
instability will exist for any value of S respectively of 
ro/^c- The change induced by choosing some different, 
more regular and differentiable core structure than our 
prescription ^) is the numerical value of the slope of the 
negative group velocity branch, within a number of or- 
der unity, and the point, as a function of ro/^c, where 
n = ric and the excitation energy becomes negative. The 
fact that around ric there are dynamical instabilities tak- 
ing place has also been recognized in jl2|, where the 
dynamical instability was investigated using the Gross- 
Pitaevskii model of a superfluid, and it was indeed found 
that ncS,c/fo = 0(1). However, what has been missed in 
this work (and others in the conventional fluid mechan- 
ical framework [^|Jll||), is that an energetical instability 
will take place right after we have crossed the dynamical 
instability region, independent from the precise value of 
the critical ric as a function of tq/^c- It is apparent from 
Fig. U that the instability will persist for any (regulariz- 
ing) model taken of the core region, as long as the energy 
density of the core is significantly less than that of the 
surrounding bulk fluid, i.e., as long as it is still signifi- 
cantly reduced compared to the bulk if we, for example, 
turn on interactions (take into account compressibility) 
inside the core. For smaller energy density differences 
between core and bulk, the value of ric will be shifted 
upwards (for a given value of ro/S.c), but the energetical 
instability will still exist. 

IV. LOW TEMPERATURE THERMODYNAMICS 
OF VORTEX OSCILLATIONS IN HELIUM II 

Up to this point, our considerations have been in 
terms of a classical vortex. Consider now the quantum 
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mechanical zero point fluctuations of a vortex line in 
the quantum fluid helium II, each with a contribution 
^hujn to the vortex (ground state) energy. If we sum 
up these contributions to the limiting value rtc, we get 
Efi = X]nL2 — 0.0035 fi^cf^"^ (cf. the area under 

the dispersion curve in Fig. H). Comparing this with the 
stationary energy Eq, see equation ( p^ ) of an undeformed 
ring in helium II, we obtain 



where, numerically. 



Ed ~ 0.014 



En 



(20) 



where the interparticle distance d — {po/m)^^/^ (^ 
in helium II) . The total quantum mechanical fluctuation 
energy in the stable modes (at zero temperature) is thus 
much less than the stationary energy Eq, of the order of 
percent. This need not be the case if we take into account 
thermal fluctuations as well. The vortex free energy may 
be written 



T{T, 6) = Eo + (3~^ ^ In (2 sinh[/3;iw„/2]) 



(21) 



n=2 



The cntropic part of the free energy, due to ring oscil- 
lations, plays an important role if the temperature is a 
significant fraction of the cyclotron energy of the core. 
We stress that the temperatures for fluctuations to be- 
come important are significantly higher if the cutoff is 
chosen well below rij,. We have also convinced ourselves 
that the absolute ratio of the oscillation free energy part 
over the stationary ring energy, \T — Eq\/ Eq, is the larger 
the smaller the radius ro is, i.e., the oscillations play an 
increasingly important thermodynamic role for smaller 
rings. 

For low temperatures, fulfilling ksT <C max [fiWn] — 
0.011 fiwc, Kelvin modes of small wavenumber and the 
modes with approximately linear dispersion around Uc 
(cf. Fig. ^, are populated |13|. In the superfluid helium 
II, where hhOc is of the order ten Kelvin, the condition 
on the temperature leads to the requirement T 100 
mK, which is feasible in experimental practice. For vor- 
tex rings of a given size and orientation, we thus expect 
two contributions to the specific heat at low tempera- 
tures, coming from the aforementioned two asymptotic 
branches of excitations on the filament. For the indi- 
cated range of mode numbers, we may approximate the 
dispersion by the Kelvin-like form 



t^/f = 7lf^^ ln[n*/n] (1 < n < nc/2) , 



where the parameters 71 and n* are 
|i I / tc ^ ^ 

71 



4ro 



ro 



Near nr. a linear law obtains: 



— 72 (j^c - n-) (ric - n <t: nc/2) , 



(22) 



(23) 



(24) 



72 ~ 0.045 WcCcAo. 



(25) 



Both of these approximate dispersion relations are valid 
for large tq/^c (small 5). The density of states for the 
w^-branch, within logarithmic accuracy, may be written 
Nk{E) ~ (4/171 ln[n*]i;)-i/2; for the w^-branch it is 
independent of the energy E, Na{E) — (^172)"^. The 
asymptotical behaviour of the vortex specific heat for low 
temperatures then assumes the form 




11.6 



(26) 



It is, as expected, proportional to the "volume", i.e. the 
circumference of the ring, and has a contribution pro- 
portional to Vt from the Kelvin-like modes and a new 
contribution proportional to T stemming from the lin- 
ear dispersion, large wavenumber branch. This last term 
gives a dependence of Cl proportional to the area 27rro^c- 
For a randomized ensemble of vortex rings, with differ- 
ent orientations and radii, we expect the indicated de- 
pendence on temperature to hold for a dilute system of 
effectively noninteracting rings. For a dense vortex tan- 
gle, coupling of the vortex rings by mutual induction will 
modify the spectrum and the above thermodynamic be- 
haviour, a problem which is left for future work. 



V. CONCLUSION 

We have derived the oscillation modes on vortex rings, 
using the canonical phase space structure of small ring 
oscillations in an incompressible, inviscid fluid and a ge- 
ometrical cutoff procedure for the core region. Beyond a 
critical wavenumber, the excitation energy becomes neg- 
ative, indicating that the vortex ring is energetically un- 
stable for perturbations on scales of short wavelengths. 
The instability relies on the energy exclusion effect of the 
deformed core, and may be interpreted in conventional 
Hamiltonian language as being due to the fact that the 
classical or quantum particle representing the excitation 
has both a kinetic energy with negative mass and a po- 
tential with negative spring constant. 

The existence and peculiar anomalous dispersion of 
propagating modes with very small wavelengths should 
have important implications for the dynamics of vortex 
reconnection events ||l^ , as well as the final stages of the 
energy cascade process in superfluid turbulence [|l^. In 
addition, we expect that scattering cross sections of the 
elementary roton excitation in helium II with vortices, 
and thus the coefficients of mutual friction between su- 
perfluid and normal components p^, will be influenced 
by the presence of low energy modes with wavenumbers 
of order the inverse core size. 
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